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Abstract 

The covariant Vlasov-Maxwell system is used to study breaking of relativistic warm plasma waves. The 
well-known theory of relativistic warm plasmas due to Katsouleas and Mori (KM) is subsumed within 
a unified geometric formulation of the 'waterbag' paradigm over spacetime. We calculate the maximum 
amplitude -Emax of non-linear longitudinal electric waves for a particular class of waterbags whose geometry 
is a simple 3-dimensional generalization (in velocity) of the 1-dimensional KM waterbag (in velocity). It is 
well known that the value of lim„^c Ema,x (with the effective temperature of the plasma electrons held fixed) 
diverges for the KM model; however, we show that a certain class of simple 3-dimensional waterbags yields 
a finite value for lim„_»c i?max, where v is the phase velocity of the wave and c is the speed of light. 

Introduction 

Considerable effort has been devoted to developing compact accelerators employing the enormous electric fields 
present in plasma wakes driven by intense lasers [1] or charged particle beams [2] (see [3,4] for recent discussions). 
Conventional accelerators operate by exciting RF microwave cavities with klystrons and use the longitudinal 
electric component of a cavity mode to accelerate bunches of charged particles for subsequent collision. However, 
it is anticipated that electric field strengths in the next generation of accelerators will be so high that the RF 
cavity walls may undergo electrical breakdown [5]. To address this issue, researchers have turned to plasma- 
based acceleration mechanisms whose field can be orders of magnitude beyond that of conventional accelerators. 
Recent years have seen the on-going development of compact sources of intense electromagnetic radiation in the 
X-ray to THz frequency range [6] that employ laser-driven plasma acceleration. Such sources promise a wide 
range of applications in medicine, material science and security. 

A sufficiently short and intense laser pulse propagating through a plasma may create a travelling longitudinal 
plasma wave whose velocity is approximately the same as the laser pulse's group velocity. However, it is not 
possible to sustain arbitrarily large electric fields; substantial numbers of plasma electrons become trapped 
in the wave and are accelerated, which dampens the wave. Indeed, the trapping phenomenon in longitudinal 
plasma waves lies at the heart of the original laser wakefield accelerator concept [1] . 

Although the evolution of a plasma wave dynamically trapping particles is complex, over the years much 
effort has been devoted to analytically understanding the upper bound ('wave-breaking limit') on the ampli- 
tude of plasma waves. Wave-breaking limits were first calculated for cold plasmas [7, 8] undergoing non-linear 
longitudinal electrostatic oscillations, and thermal effects were later included in non-relativistic [9] and rela- 
tivistic [10-12] contexts. The results for the cold plasma are uncontroversial, but recent discussion [13-15] has 
uncovered difficulties with establishing an agreed analytical description of longitudinal wave-breaking in warm 
plasmas; in particular, it has been noted that different plasma models based on different assumptions yield 
different results. Models of non-linear plasma waves near breaking are approaching the limits of their domain of 
applicability, and different models exhibit different wave-breaking limits. Although recent experiments [16-18] 
operate in the 3-dimensional 'bubble' (or 'blow-out') regime [19] and exploit transverse wave-breaking [20], 
recent work [13-15] has rekindled interest in the theory of longitudinal wave-breaking. 

Recent discussion [13-15] includes comparison of the behaviour of the relativistic 'waterbag' model [10,21] 
due to Katsouleas and Mori (abbreviated as KM) and a warm plasma model [12] due to Schrocder, Esarey and 
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Shadwick (abbreviated as SES) employing velocity moments of the 1-particle plasma electron distribution. The 
KM and SES models yield different results for the maximum amplitude of non-linear electrostatic oscillations in 
the limit u — > c with the temperature of the plasma held fixed {v is the phase velocity of the plasma wave with 
respect to the laboratory frame) . The KM maximum electric field diverges logarithmically in 7 = 1/ a/1 — jc? 
as w ^ c, whereas the SES maximum electric field tends to a finite value as w ^ c (with the initial plasma 
temperature held fixed in the limit v c). Employing velocity moments of the 1-particle plasma electron 
distribution, SES require that the distribution remains narrow in velocity spread whereas the KM approach 
employs a particular waterbag solution to the Vlasov equation. The disagreement of the two approaches has 
been attributed to the waterbag's piecewise constant structure and lack of a tail [14]. 

The KM model is formulated over 2-dimensional spacetime and SES employ a line distribution in longitudinal 
velocity to simplify their field equations on 4-dimensional spacetime. Neither model admits a plasma electron 
distribution with a non-vanishing transverse thermal spread. Thus, a theory of waterbags over 4-dimensional 
spacetime was recently developed [22, 23] to permit analytical investigation of wave-breaking as a function of 
the 3-dimensional shape (in velocity) of the plasma electron distribution. In the following we cast the KM field 
equations in a form comparable with those of the waterbag over 4-dimensional spacetime and, for the first time, 
give a unified presentation of the derivation of wave-breaking limits in the KM model and our waterbag model. 
We conclude with a comparison of the predictions of a particular class of our waterbags, the KM model and 
SES model. We find that the results of our present approach have more in common with the SES model than 
the KM model. 



1 Vlasov-Maxwell equations 

The brief summary of the Vlasov-Maxwell equations given below establishes our conventions; further details 
may be found in [23,24]. We employ the Einstein summation convention throughout this article. Latin indices 
a, 6, c run over 0, 1, 2, 3 and units are used in which the speed of light c = 1 and the permittivity of the vacuum 
eo = 1- 

Let (x") be an inertial coordinate system on Minkowski spacetime [M^g) where is the proper time of 
observers at fixed Cartesian coordinates (x^, a;^, x'^) in the laboratory. The metric tensor g has the form 

g = ijab dx" dx'' (1) 

where 

-lifa = 6 = 0, 

Vab= { 1 if a = 6 / 0, (2) 
if a / 5. 

Let (x", x'') be an induced coordinate system on the total space TM of the tangent bundle {TM,Il, M) and 
in the following, where convenient, we will write x instead of x"" and x instead of i''. For notational simplicity, 
we will not distinguish between a point in a manifold and its coordinate representation. 

The total space £ of the sub-bundle {£,Il\£,M) of {TA4,Il, M) is the set of timelikc, future-directed, unit 
normalized tangent vectors on A4, 

£ {{x,x) eTM\(p = a.nd x° > 0} (3) 

where 

V? = 77ati°i''-f 1. (4) 

Plasma electrons are described statistically by a 1-particlc distribution / on TA4 which induces a number 
4-current vector field N, 

N-{x)= [ xyLx#l, (6) 
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where Ex = (njf) ^{x) is the fibre of {£,I1\£,M) over x e M. The 3-form on TM is induced from the 

4-form #1, 

#1 = di" A di^ A di^ A c^i^ (7) 
and the dilation vector field X, 

on TM, where lx is the interior product on forms. It may be shown 

^ ^ ^ dx'^ A dx^ A dx^ (9) 

where \x\'^ = [x^Y + (i^)'^ + (i^)^ and ~ denotes equality under restriction to £ by pull-back. The above are 
specialised to inertial coordinates on Minkowski spacetime; their form in a general coordinate system may 
be found in [23,24]. 

We are interested in the evolution of a plasma over timescales during which the motion of the ions is 
negligible in comparison with the motion of the electrons. We assume that the ions are at rest and distributed 
homogeneously in the laboratory frame. Their worldlines are trajectories of the vector field A'ion = nio^d/dx^ 
on M. where rijon is the ion number density (a positive definite constant) in the laboratory frame. The Maxwell 
equations are 

dF = 0, (10) 
d*F = -qi^N + qirN^n (11) 

where F = ^Fat dx°' A dx'^ is the electromagnetic 2- form and q is the charge on the electron (with q <Q). The 
Hodge map * is induced from ([T]) and the volume 4-form ^1, 

*1 = dx° A dx^ A dx^ A dx^ , (12) 

on M.. The 1-forms N, TVion are the metric duals of the vector fields N, iVion respectively, i.e. the 1-form Y 
satisfies Y{Z) = g(Y, Z) for all vector fields Z on A^. 

The scalar field / satisfies the Vlasov equation, which may be written 

'n^S;')-" (13) 



,dx°- m dx''^ 

on £, where F'^Y is the vertical lift of F\ ^ r]'"'Fca from M to TM, 

F^Y{x,x) = F\{x), (14) 

and m is the electron rest mass. 

The equations of motion for a waterbag distribution are readily motivated via a global expression of the 
local Vlasov equation p3|l . Introduce the Liouville vector field L, 



on TM and the 6-form lj, 

cj = tiiA-(*l'^A#l) (16) 
where the 4-form *1^, 

= dx° A dx^ A dx^ A dx^ , (17) 
is the vertical lift of the spacetime volume 4-form Tk-l from M to TM . It can be shown 

dujc^Q (18) 
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and the Vlasov equation (|T3)) can be written 

difLu) ~ 0. (19) 
Thus, it follows 

/ d(/c^)=0 (20) 
where S is a 7-dimensional region in £ and using Stokes' theorem on forms (see. for example, [25,26]) we obtain 
fio^O (21) 



I OB 

where dB is the boundary of B. 

1.1 Waterbag distributions 

We consider distributions for which / = a is a positive constant inside a 7-dimensional region hi d £ and 
/ = outside. In particular, we consider U to be the union over each point x ^ Ai oi a, domain Wx whose 
boundary dWx in £ is topologically equivalent to the 2-sphere. Such piecewise constant distributions are called 
'waterbags'. 

Choosing B in pT|) to be a small 7-dimensional 'pill-box' that intersects dWx and evaluating the integral in 
the limit as the 'height' of B tends to zero, we recover a jump condition on Jlo that leads to 

dXhujc^Q atA = (22) 

where A = is the union over x of the boundaries dyVx- If A = is the image of the embedding map S, 

Y.:My.S^ £ 

(x,0 ^ {x,x = V^[x)), (23) 
where ^ € S'^ has coordinates (^\^^), then it follows [23] from (0 [El [TH]) that ([22]) is equivalent to 

{yy^%- S^iy^F) (24) 



Here, and Vl^ are families of vector fields and 2-forms on M. respectively, where 
dx" 



= y^iL' (25) 



= -g^dXa A ^dxt. (26) 

with dxa = riabdx^ . Note that since the image of E lies in £ it follows that, for each ^ G 5^, is timelike, unit 
normalized and future-directed: 

g{V^,V^) = -l, g{V^,^)<0. (27) 

We adopt ([M)) as the equation of motion for the waterbag boundary dWx- 
It may be shown that a particular class of solutions to ((M)) satisfies 

Ttl — 

F=-dV, (28) 

and using (jlip we obtain the field equation 

di,dV^ = ~-^{i^N ^ *iVio„) (29) 

on M with the condition that dV^ is independent of ^. For simplicity, we have neglected the direct contribution 
of the driver (laser pulse or particle bunch) to the total electromagnetic field in 
For a discussion of solutions to that do not satisfy ([25]) sec [23] . 
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M2 Pi Pll 



Figure 1: Illustration of a waterbag distribution over 2-dimGnsional spacctimc A^2 with pi ~ (ti,zi) and 
Pll = {hi, zii). The shaded region is the interior of the waterbag (where / is non-zero), and the 2-vector fields 
{V+, V-} determine the boundary of the waterbag. 



2 Electrostatic oscillations on 2-dimensional spacetime 

Before analysing ([?7l UHl) further it is useful to briefly discuss their analogue on 2-diinensional spacetime for 
facilitating comparison with the approach adopted in [10,21]. Although formulated on 4-dimensional spacetime, 
equations (|271 129p have a similar structure for any number of dimensions. In particular, we now consider 2- 
dimensional Minkowski spacetime (A^2,ff) 

g = -dt (E)dt + dz(E) dz, (30) 
-kl^dt Kdz (31) 

where (t, z is a Cartesian coordinate system in the laboratory inertial frame. An induced coordinate system 
on TM.2 is {t, z, t, z) and the 2-form #1 and dilation vector field X over TM.2 are 

#1 = A dz, (32) 

^ = 4 + (33) 
ot oz 

Furthermore, ^ is now an element of the 0-sphere {+, — } and = 1 is a constant 0-form. Thus, the analogue 
to dMl) is 

Vv^V+-^iv^F^Q, (34) 
m 

Vv^Vl - —Lv^F ^0, (35) 
m 

where satisfy the conditions 

g{V+,V+)^-l, g(F+,^)<0, (36) 
g{V-,V-) = -l, g(F_,^)<0, (37) 

and the only non-trivial Maxwell equation for the 2-form F is 

di.F = -qi.N + qi.N^n (38) 



^We use (t, z) rather than (x") to distinguish coordinates on 2- and 4-dimensional spacetimes. 
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where TVjon = riiond/dt is the ion number 2-current and F = Edt A dz where E is the electric field along the 
z-axis. 

On the unit hyperbola bundle £2, t ~ vl + i^ and the components of the electron number 2-currcnt 
N — N^d/dt + N^djdz corresponding to (O are 



= f f{t,z,i,z)dz ^ a(Y+ -Y^^, 



N'= I ^==f{t,z,t,z)dz^a{^/l + Yf-Jl + Y^), (39) 

where 



|0, i < y„ or i > y+ 

with a a positive non-zero constant and 1^-} 0-forms on A^2- The 2-velocity fields V-} satisfy 
and it follows 

TV = a*(V^ - V^). (42) 
See figure [TJ 

Unlike their 4-dimensional analogue, which may include transverse electromagnetic fields, p4l35p are uniquelij^ 
solved by 

dV± = —F (43) 
m 

and using ([55]) it follows 

_ „2 ^ 

di,dV± = - — (i.N - ★TVion ) (44) 
m 

subject to the condition dV+ = dV-. 

Alternatively, one may follow the approach adopted in [10] by casting the above as a warm fiuid. The type 
(0,2) stress-energy-momentum tensor Tfluid of the electron fluid is 

/ /■ z^i^f \ d d 



where Greek indices ji, v run over 0, 1 and z^ ^ t, z^ = z, i° = VT+i^, z^ = z. It can be shown that Tfluid 
induced by the above waterbag distribution can be expressed entirely in terms of the proper number density n 
of the electron fluid, the electron fluid's bulk 2-velocity U and the spacetime metric: 

Ifiuid = {p + p)U (E) U + pg (46) 

where g is the inverse metric tensor 

d d d d 

at ot oz az 

and 

U^^2==Z, Z = \{V+ + V.), (48) 
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N = nU, n= ^J-g{N, N), (49) 



^Proper incorporation of transverse fields requires at least 2 spatial dimensions. 
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with the equations of state 



P 



ma 




sinh-i^ 
2a 



(50) 
(51) 



The equation of motion of the electron fluid, 

ip + p)VuU ^qmuF - LuidpAU), (52) 
follows from the zero divergence of the sum of Tfiuid and the Maxwell stress-energy-momentum tensor where 



g{U,U) = -l, 



.(c/,|)<o. 



(53) 



It should be stressed that the warm fluid model ([501 [5T1 [5^ is equivalent to [551 [Ml [57)) . Thus, 
I35[ 1361 I37p may be replaced by an equivalent field theory expressed in terms of a finite set of moments of 
/ on 2-dimensional spacetime. However, the situation is more complicated for waterbags over 4-dimensional 
spacetime where second, and higher, order moments of / in i; are not, in general, easily expressible in terms of 
zeroth and first order moments of /. In general, the moment hierarchy does not automatically close. 

We will now obtain a non-linear ordinary differential equation describing 1-dimensional electrostatic oscilla- 
tions and determine an expression for the wave-breaking limit of this model. Derivation of wave-breaking limits 
starting from Tfluid and the equations of state ([501 [STjl may be found in [10,21]. However, we will work directly 
with ([501 [57[ [m mi) to facilitate comparison with our model on 4-dimensional spacetime. 

Let all field components with respect to the laboratory frame (dt, dz) be functions oi Q ~ z — vt only (the 
'quasi-static assumption'), where < w < 1, and let (e^,e^) be the basis 



vdz — dt. 



dz — vdt. 



(54) 



The coframe (76^,76^) is an orthonormal basis adapted to observers moving at velocity v along z (i.e observers 
in the 'wave frame') where 7 = (1 — u^)~^/^ is the Lorentz factor of such observers relative to the laboratory. 
So, 7e^(Afion) = —jnionV is the component of the ion number 1-current in the wave frame. 
In the basis (e^, e^), V± can be decomposed as 



V± = (M(C) + ^±)ei+^±(C)e' 



(55) 



where {A+, A-} are constant. Note that this is the most general decomposition compatible with equation ([^5)) 
and the quasi-static assumption. 
Solving ([Ml 1311) for V'i gives 

'/'± = (m + A±)'-7' (56) 

and additional physical information is needed to fix the sign of ^p± . Here, we demand that all electrons described 
by the waterbag are travelling slower than the wave so ip± 

1/2 



V(m + ^±)^ -7^ and ^ is 



V± ^{f^ + A±)e' - {^{fi + A±) 

Substituting (j55p into equation 
^ 1 m d/i 

72 5 dc' 



7 

yields 



and equation 
1 d'^p 

i^~d^ 

with the algebraic constraint 

„2 

A+-A^^ 



yields the nonlinear oscillator equation 
2 



<f 2 

7 "ion 

m 



ii7 " 



< 0. 



v/(m + ^+)2-72-v/(m + ^-)'-7^ 



a 



(57) 



(58) 



(59) 



(60) 
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2.1 Electrostatic wave-breaking 

In the wave frame the relativistic energies of the two ends of the waterbag are m(/i + A_|_)/7 and m(/i + A-)/^ 
respectively, and since m{p. + A+)/7 > m it follows ^ + A± > 7. Using (pU]) . fi + A- > /i + A+ and hence 
H + A^ > 7 implies /i + A_ > 7. Thus, /i + A± > 7 may be reduced to ^ > //wb where 



(61) 



Alternatively, one may arrive at the same conclusion by inspecting the right-hand side of (|59p and using fJ.+A± > 
(which follows because V+ and V- arc future-pointing). Thus, there is an upper bound on the amplitude of 
oscillatory solutions to (|59p . which leads to an upper bound E„iax on the electric field E (the 'wave-breaking 
limit' of this model). 

During an oscillation E vanishes when dfi/d(^ vanishes and is at a maximum when \dfi/d(\ is at a 
maximum (see (j58p and note q < 0). A maximum of \dfi/d(^\ occurs at values ^0 of C where ^(Co) equals the 
oscillator equilibrium ^cq- Furthermore, for the maximum amplitude oscillation d^/dC, vanishes when = /^wb- 
An upper bound i^max on the magnitude \E\ of the electric field is obtained by evaluating the first integral of 
([5^ between /i = /i^b and /i = /icq. 

Without loss of generality, we can choose the split between /x and A± such that 



A+ = -A^ = -a, 
Using (HHUnH) it follows 



2a 



Mcq + Mwb + 



1 V 

2a 



where figq is the equilibrium solution to ((5^ . which satisfies 
2a 

V 



(/icq + a)2 - 7^ - V (Meq - «) 



72 d/i 



(62) 



(63) 



(64) 



The constant a is fixed in terms of an effective temperature Teq associated with the oscillator equilibrium 
Heq- Noting that n = njon when the waterbag is in its equilibrium state (/it — /^eq), and assuming riion ^ 2a in 
(fCT]) it follows 



Pcq 



3(2a)2' 
Introducing Toq via poq 



nionksTcq, where ks is Boltzmann's constant, we find 



^ion 



oq 



7^u 



3fcsTc 



oq 



(65) 



(66) 



where ([5^ has been used. Hence, nion «C 2a means that the thermal energy of the electron fiuid in the oscillator 
equilibrium state is much less than the rest mass-energy of the electron. 

The wave-breaking limit -Emax can be readily analysed for 7 3> 1 via asymptotic approximation in a small 
parameter e, 



_ 7 _ 2a 1 

a TT-ionf 7 ' 

where (j62|) has been used. Employing (j6H [63l [64|) it follows 

2 r 1 

^max = 2mniona - Acq + Awb 



1 

2'' 



[Acq + 1]^ 

Awb = 1 + e 



[Acq - 1]' 



1, 



(67) 



(68) 

(69) 
(70) 
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where, /t = /i/a. To proceed further, we express v in (j68l I69|) as a function of e and a parameter b that 
characterizes the effective temperature of the oscihator equihbrium distribution. Using ([62l [67]) it foUows 



^ (71) 



where b is 

'-If^ (72) 

The dominant e dependence (as £ — > with & held fixed) of iJ^ax arises from the second term in the integrand 
in (j68p and may be extracted by expanding the integrand with respect to e and integrating each summand over 
jj,. Since, for > e > 0, 

n— 2 

where Cn are constants, and inspection of (j69p reveals 

Acq = /i(e') = MO) + /i'(0)£2 + 0(e4) (e^O), (74) 
we find, for h{0) > 1, 



// 1 1 \ '^"'^ 1 

V / ,1—9 



-£2 



+ -e2hr(£)+0(£2) (e^o) (75) 

Awb 



where (|70p has been used. Furthermore, it foUows from (|69p that an asymptotic approximation for /i(0) in smaU 
b leads to 

M0) = i+O(l) (6^0). (76) 
Thus, IMl) yields 

:^«-i™nione'ln(£) (77) 
a L 

for e,b <^ 1. Introducing the effective temperature Toq using (|621 1661 1^7]) and noting v ~ 1, we obtain 



^ax^:^— r^J:Tl-^ln 7a/-^ , (78) 



2 g2 W 3fc^r, 



oq 



where ujp = \/ Uionq^ / [meo) is the plasma frequency and the speed of light c and permittivity of the vacuum £0 
have been restored. Equation (j78p was obtained as a lower bound on Ef^-^^^ in [13]. 



3 Longitudinal electrostatic oscillations on 4-dimensional spacetime 

We now consider longitudinal electrostatic waves on 4-dimcnsional spacetime by closely following the above 
description on 2-dimensional spacetime. 

As before, we adopt the 'quasi-static assumption'. We seek a waterbag Wx axisymmetric about whose 
pointwise dependence in Minkowski spacetime M is on the wave's phase ( = — vx'^ only, where < w < 1. 
As before, the following results are applicable only if the longitudinal component of in the wave frame is 
negative (no electron described by Wx is moving faster along x^ than the wave) . 

Decompose in the wave frame as 

Vi = HO + M^)] + ^i^\0 + i?sin(^i) cos{^^)dx^ + Rsm{^^) sm{^'')dx^ (79) 
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for < ^1 < TT, < < 27r where i? > is constant and 

vdx^ ~ dx° , e^^dx^-vdx^. (80) 



Here, (76^, 76^, dx^ , dx'^) is an orthonormal basis adapted to the wave frame, with 7 = l/\/l — f^. In the wave 
frame the relativistic energy of = mV^ is m(/i + A(^i))/7 and it follows that /i + ^(C^) > 0. Furthermore, 
using dlTl [711) it follows 



^ = -^/[/i + - 72[1 + i?2 sin2(ei)], (81) 

where the negative square root is chosen because no electron is moving faster along x^ than the wave, and we 
obtain 



^^>-A{e) + l^Jl + R'sin\e)■ (82) 
Substituting ([TQ]) into equation ((28l) leads to 

^ m da r, , , ^ 

F=--^e^Ae\ 83 

and dHISaillllHIl) yield 

^ ~^nion7'--27ri?2a / f[;, + A(^i)]2 - ^^[^ ^ ^2 ^^^^2(^1)^ ^ sin(ei) cos(^i) (84) 

m TO J \ / 
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(c.f. equation ([59]) ) and 

TT 

27ri?2 / ^(^1) sin(^i) cos(Ci) = (85) 


(c.f. equation ([60]) ) where a > is the value of / inside Wx- 

The form of the 2nd order autonomous non-linear ordinary differential equation (j84p for /i is fixed by 
specifying the generator A{£^^) oi dWx subject to the normalization condition ([55)1 . 

3.1 Electrostatic wave-breaking 

The form of the integrand in ensures that the magnitude of oscillatory solutions to (|84p cannot be arbitrarily 
large. For our model, the wave-breaking value /^wb is the largest fi for which the argument of the square root 
in ([84]) vanishes. 



Aiwb = maxj - A{^^) + j^Jl + R^ sm\e) < < 7r| 



(86) 



because ^ < fi^rh yields an imaginary integrand in (|84p for some 

The electric field has only one non-zero component E (in the x^ direction) . Using F ^ E dx^ A dx^ and ([80 
[551) it follows 



E^--,^ (87) 



TO 1 d/i 

q 72 dC 

and the wave-breaking limit i^max is obtained by evaluating the first integral of (|84p between /i^b where E 
vanishes and the oscillator equilibrium /igq of /i where is at a maximum. Using (|85p to eliminate a it follows 
that /ieq satisfies 

i 1 ^(^1) sin(Ci) cos(ei) d^i = j (^[/icq + A(Ci)]2 - 7'[1 + i?' sin2(Ci)]^ s\n{e) cos(^i)dC 
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with 



TT 



(89) 



since a,v > 0. Equation (|84p yields the maximum value i?niax of E, 



El 



2mn;, 



Mcq + Mwb+ — 



[/i + - 7'[1 + i?' sin2(ei)] sin(ei) cos(^i)dei d/i 



1/2 



(90) 



To proceed further we need to choose the generator A{S}) of the waterbag distribution. It turns out that 
even the simple choice 



^(e^) = -acos(ei) 



(91) 



for A{^^), where a is a positive constant, leads to a wave-breaking limit -Emax with interesting behaviour, as we 
now show. 

Using dnni) it follows 



2771711, 



iv 
2a 



Mcq 1 



1/2 



(92) 



where x = — cos(^-^), equation ([55]) yields 



\l j {[^,^^ + axr~l'[l + R\l~x')]j'\dx^l 



2a 

and equation (|86p may be written 



-ax + l^l + R^{l-X^) 



-1 <X< 1 



Examination of reveals that two classes of waterbag arise according to whether or not the function 
X^-ax + l^/l + R^{l-X^) 



(93) 



(94) 



(95) 



has a turning point in the interval [—1, 1]. Examples of the two classes are shown in figures [5] and [21 In each 
case, the plasma wave breaks when the uppermost part of the distribution achieves the phase velocity of the 
plasma wave (i.e. the longitudinal component tl) of in the wave frame vanishes). Wave-breaking limits for 
the class in figure [5] have been calculated previously [22,23] and here we focus on waterbags of the type shown 
in figure [3l 

3.1.1 Calculation of the maximum electric field 

The parameters a, i?, 7 are chosen to satisfy 



1 + i?2 



R \ a? + 72i?2 



> 1, 



(96) 



ensuring that (j95|) does not have a turning point in the interval [—1, 1]. Hence, the wave breaks when the tip 
X — cos(O) = — 1 of the waterbag achieves the phase velocity of the plasma wave. Using ([M]) it follows 



Mwb = a + 7 



(97) 



11 



x^ 




Figure 2: Two illustrations of the 4- velocity dependence of a particular 'bowl' waterbag. The axis of symmetry 
is aligned along x^. The maximum electric field amplitude is achieved during the oscillation in which the top 
of the waterbag (a circle) grazes the phase speed of the wave. 




• 1 -2 



Figure 3: An illustration of the 4-velocity dependence of a particular 'gourd' waterbag. The axis of symmetry 
is aligned along . The maximum electric field amplitude is achieved during the oscillation in which the tip of 
the waterbag grazes the phase speed of the wave. 
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which is formally identical to the wave-breaking limit of /i for the waterbag over 2-dimensional spacetime. This 
is quite different from the value of /^wb for waterbags of the type shown in figure [2] (see [22,23]). 

Following a similar method to that used in section [21 we now evaluate (I92p for 7 3> 1. Introducing /i = /i/a 
in dH m nil) leads to 



-E'max = Smnionfl 



Mcq 1 



1/2 



Xdp-dx 



1/2 



= 1, 



Awb — 1 + e, 

where, using (|85l) . 

_ 3nion7'^t' 
47ri?2Q! 

e = — 



a SriionT^' 

Thus, it follows e as w 1 and we determine an asymptotic approximation for i?max in e as e — > 0. 
Expansion in e of the integrand in ([551) yields 



(98) 

(99) 
(100) 

(101) 
(102) 



1/2 



A + x- 



2(A + x) 



)i=2 



(A + x) 



where the c„ are numerical constants. Using (|103p . the integral over A in ([Ml) leads to a summand proportional 
to 



fn 



2\\n 1 



(l + fl^(l-X^))" {l + R'jl-x')) 

(Aoq+X)'"-' (l+£ + X)2"-' 



n > 2 



(104) 



where (jlOOp has been used. 

Inspection of (|99p suggests an approximation for Acq(£) of the form 

Acq(£) = /i(e') - HO) + h'{0)e^ + O(e^) (e ^ 0). 

Using ([TU^l) it follows 

1 3721011 a 



and ([Ml) leads to 

3 1 + R^l~x^) 

Acq(O) + X 



Xdx = b 



Thus, Aoq(O) niay be approximated as 



1 / 2/72 

Acq(o) = -^i + — + o(i) ib^o). 

Repeated integration by parts in p04p leads to 
/„=0(£3-2«) (e^O), n>2 



(105) 

(106) 
(107) 

(108) 

(109) 
(110) 
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and we obtain the asymptotic approximation 



1 Acq . -^^2 



[A + x]'-e1l + i?'(l-x")] Xdxdfi 



-1 l+£ 



1 

^(Aeq - Awb) - J[l + R'il x')] In ]xdx + 0{e') {e ^ 0) 



kAeq - Awb) - J[i + R'ii - x')] In (^TTT^)^ + ^'^"^ 

(111) 



Thus 



-1 



fe^(l-Acq(0)) + - J[l + R'{l-x')]lni^-^-^^^]xdx}+OieHns) (e 0) (112) 



and retaining lowest order terms in e. 6, i? yields 

1 

^Lx ~ ^mnionae^ /[I + i?'(l - x')] ln(l + x)xt^X, 



-1 



« -mriion^ (113) 

where (|107p has been used to eliminate aR^. 

Numerical validity of the above approximation is supported by figured! The solid curves are obtained by 
numerically integrating ()98m00p and the dashed lines are obtained using pi3p . It is clear that pi3p yields a 
good approximation to Smax for large 7. 

In order to compare (|113p to expressions for i?niax obtained elsewhere [10-12], it is useful to express (|113p 
as a function of effective temperature. The electron proper number density n = riion when ^ = and we 
eliminate b in favour of an effective longitudinal temperature Tygq defined as 

I^B "■ion 

where ks is Boltzmann's constant and p\\cq is the longitudinal pressure associated with the oscillator equilibrium 
ji = /Zoq. The longitudinal pressure is p\\cq ~ %?q where the stress-cnergy-momentum tensor 7^q has components 

T^f = ma [ x^xhx^l (115) 

with Weq the support of the waterbag distribution jj — /icq (the choice of fibre is unimportant as the distribution 
associated with /icq is independent of C)- 
Since 

x^^ dx^ A dx^ A dx'^ . , , 1 ^9 . , _i /i^ M , . 1 ,.9! , 

dx A dx > (116) 



where i^^ = (i^)^, x^^ = (i^)^, i^^ = (i^)^ and /? = Vl + + i^^, using (|115p and Stokes' theorem on 
forms, it follows 



iiV/3'+i^'- sinh-i (^] dx^ A dx''. (117) 
2 2 V ' 
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Using (|791 [50]) , components of the oscillator equilibrium waterbag = V^^^ are 



i° - A^eq - a cos(ei) - t- \/ [Meq - a 003(^1 )] ^ - 7^ [1 + sin^ (^i )] , (118) 
= i?sin(^i)cos(^^), (119) 
= Rsm{^^)sm{S,^), (120) 



±3 = v[^l,^ -acosie)] - VlA^cq - acos(ei)]2 - ^^[1 + sin^C)] (121) 



and it follows 

^ = i;[Acq - cos(Ci)] - ^[Acq-cos(^i)]2-£2[l + i?2sin2(ei)] 

- i52[Acq(0) - cos(ei)] + ^ (122) 

2 2[^eq(0) - COs(Cl)]y 

to lowest order in e and i?. Using p091 1122p it follows 

— ^eH''cos{e) (123) 
a 

to lowest order in e, 6 and i?. Furthermore, (|102U107]) yield ae^6^ = b/v b and so 

i:3«6cos(^i) (124) 
to lowest order in e, 6 and R. Hence 

^^/32 + i32_ 1^2ginh-i (^^^ Ib^cos'ie) (125) 
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to lowest order in e, b and R and (|117p yields 



Plleq 
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(126) 



Equations (fTTl fTTH [T^ yield 




( 



) 



1/2 



max 



20/cBT||eq 



e,5,i? < 1 



(127) 



where Wp = \/nionq'^ /{meo) is the plasma frequency and the speed of light c and permittivity of the vacuum Eq 
have been restored. 

Conclusion 

Equations ([751 I127P indicate that waterbags over 2-dimensional spacetime and 4-dimensional spacetime can 
behave quite differently. Equation (|127p is independent of 7 but ((75|) diverges as 7 ^ 00, and this difference 
in behaviour arises because the logarithmic singularity in the integrand in (|113p is integrable. Moreover, the 
^||eq^^ behaviour of the asymptotic form of i^max for ksT^^cq ^ is very similar to the results of SES [12] 
and others [11] in the limit tj — > c. 

Direct comparison of our results and those of SES follows by setting the transverse vector potential A± to 
zero in the SES model, thereby neglecting the overlap of the electromagnetic field of the driver (laser pulse or 
particle bunch) and the wave. The approach followed by SES begins with covariant field equations, induced from 
the Vlasov equation, that couple the zeroth, first and second order centred moments (in i°) of the 1-particle 
distribution / with the electromagnetic field. SES then assume that the 1-particle distribution / (restricted by 
pull-back to the unit hyperboloid) may be approximated f ~ h{x'^ , , x'^)S{x^)d{±'^) where S is the Dirac 
delta function. A covariant measure of the total thermal spread is given by the magnitude of the ratio of the 
trace of the second order centred moment and the zeroth moment. SES assume that the third order centred 
moment is 0{e^) and can be neglected relative to lower order moments. 

One could develop a similar argument to that given by SES based on moments of a prescribed 3-dimensional 
waterbag with narrow velocity spread, rather than the line distribution employed by SES. However, nuances in 
the shape of the waterbag would be lost; for example, we would not know that merely the tip of the waterbag 
grazes the wave's phase velocity (see figure ^ during the maximum amplitude oscillation. This could be 
important because, as noted earlier, longitudinal wave-breaking is associated with the trapping of considerable 
numbers of particles in the wave (see [15] for a discussion), and our present model neglects trapped particles. 
Thus, we expect that £^max calculated here is a lower bound on the maximum electric field obtained when 
trapping is accounted for. 

In conclusion, we have shown that it is possible to construct 3-dimensional waterbag distributions that lead 
to a maximum electric field amplitude whose asymptotic behaviour is similar to that of the SES model as w — > c 
(with effective temperature held fixed in the limit v c). 

Acknowledgements 

We thank RMGM Trines for useful discussions. We acknowledge EPSRC for financial support. 

References 

[1] T Tajima and JM Dawson, Phys. Rev. Lett. 43 (1979) 267 
[2] P Chen, et al, Phys. Rev. Lett. 54 (1985) 693 
[3] V Malka, et al, Nat. Phys. 4 (2008) 447 
^We have changed the notation used by SES to avoid conflict with our own. 



16 



[4] A Caldwell, et al, Nat. Phys. 5 (2009) 363 

[5] W Wuensch, Proc. EPAC 2002 134 

[6] HP Schlenvoigt, et al, Nat. Phys. 4 (2008) 133 

[7] Al Akhiezer and RV Polovin, Sov. Phys. JETP 3 (1956) 696 

[8] JM Dawson, Phys. Rev. 113 (1959) 383 

[9] TP Coffey, Phys. Fluids 14 (1971) 1402 
[10] T Katsouleas and WB Mori, Phys. Rev. Lett. 61 (1988) 90 
[11] JB Rosenzweig, Phys. Rev. A 38 (1988) 3634 

[12] CB Schroeder, E Esarey and BA Shadwick, Phys. Rev. E 72 (2005) 055401 

[13] RMGM Trines and PA Norreys, Phys. Plasmas 13 (2006) 123102 

[14] CB Schroeder, E Esarey and BA Shadwick, Phys. Plasmas 14 (2007) 084701 

[15] RMGM Trines and PA Norreys, Phys. Plasmas 14 (2007) 084702 

[16] SPD Mangles, et al. Nature 431 (2004) 535-8 

[17] CGR Geddes, et al. Nature 431 (2004) 538-41 

[18] J Faure, et al. Nature 431 (2004) 541-4 

[19] W Lu, et al, Phys. Plas. 13 (2006), 056709 

[20] T Esirkepov, et al, Phys. Rev. Lett. 96 (2006) 014803 

[21] WB Mori and T Katsouleas, Phys. Scr. T30 (1990) 127 

[22] DA Burton and A Noble, AIP Conf. Proc. 1086 (2009) 252 

[23] DA Burton, A Noble and H Wen, II Nuovo Cim. C 32 1 (2009) 1 

[24] J Ehlers in General Relativity and Cosmology, Proceedings of the International School of Physics "Enrico Fermi" 
47, (Academic Press, New York and London, 1971) 1 

[25] DA Burton, Theoret. Appl. Mech. 30 (2003) 85 

[26] IM Benn and RW Tucker, An Introduction to Spinors and Geometry with Applications in Physics (Adam Hilger, 
Bristol and New York, 1987) 



17 



